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Abstract 

We present calculations of the spin-averaged pair distribution function g(r) in a homogeneous gas 
of electrons moving in dimensionality D = 3 or Z) = 2 at finite temperature. The model involves 
the solution of a two-electron scattering problem via an effective potential which embodies many- 
body effects through a self-consistent Hartree approximation, leading to two-body wave functions 
to be averaged over a temperature-dependent distribution of relative momentum for electron pairs. 
We report illustrative numerical results for g{r) in an intermediate-coupling regime and interpret 
them in terms of changes of short-range order with increasing temperature. 
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I. INTRODUCTION 



Many of the electron-electron interaction effects in condensed-matter and astrophysical 
plasmas can be understood by reference to the homogeneous electron-gas model. A central 
role is played by the electron-pair distribution function g{r), which describes how short-range 
exchange and correlations enter to determine the probability of finding two electrons at a 
relative distance r. Most of the theoretical and computational studies of this model have 
been concerned either with the limit of complete quantum degeneracy at zero temperature 
or with the classical limit at very high temperature Q]- In a number of physical situations, 
however, the state of the electrons is in neither of these limiting regimes. Some studies 
of the equation of state of the three-dimensional {3D) electron gas in the temperature- 
density plane were reported most notably by Ichimaru and coworkers and by Macke and 
coworkers P|. In fact, the excess free energy of the electron gas at any temperature T can 
be calculated from the temperature and density dependence of g{r) by integration over the 
coupling strength parameter. 

Accurate and theoretically well founded methods have been developed in recent years 
for the calculation of q(r) in both the 3D and the 2D electron gas at T = 0. Here we 

n 

are specifically concerned with the approach originally proposed by Overhauser j5|] for the 
evaluation of the value g{0) of the pair function at contact. Writing g{r) in the paramagnetic 
state as the average of the distributions for parallel- and antiparallel-spin electron pairs, 
g{r) = + (7||(r)]/2, one has that g^^{0) vanishes on account of the Pauli principle and 

g^l{0) is solely determined by two-body scattering events. In Overhauser's model g^i^O) was 
obtained from the solution of a Schrodinger equation for the relative motion of two electrons 
with antiparallel spins in an effective interaction potential which also accounts for many- 
body effects from the electron gas. This approach was later extended to evaluate g{r) at 
finite r as a sum of contributions from two- body wave functions ( "geminals" ) , the geminal 
weights being determined by the distribution of relative momentum for a pair of electrons 
in the electron ga.0flaa. 

In the present work we extend to finite temperature the approach of Davoudi et al. [7(, 
which was based on a self-consistent Hartree model for the effective scattering potential. 
Section 2 summarizes the method used for the evaluation of the pair wave functions from 
the solution of an effective two-body Schrodinger equation in dimensionality D = 3 and 
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D = 2, while Section 3 presents the calculation of the temperature-dependent geminal 
weights needed to construct the pair distribution function. Some illustrative numerical 
results for g{r) are presented and discussed in Section 4, and finally Section 5 offers a brief 
summary and some perspectives for future work. 

II. TWO-BODY WAVE FUNCTIONS 

The thermodynamic variables that characterize the states of the gas in dimensionality D 
are the electron density n and the temperature T. It is convenient to introduce the Wigner- 
Seitz density parameter by the usual relation r^aB = [D/ (ti^Id)]^^^ , where is the Bohr 
radius and fi^ is the solid angle in D dimensions {Q^ = Att and ^2 = Stt). We shall measure 
the thermal energy in eV units by introducing a reduced temperature B = kBT/eV. 

As in the gas at zero temperature , we aim to solve the two-electron scattering problem 
with some effective interaction potential V{r) in order to determine the wave functions 
\I^j^^j^(r, R) and \E'j^^j^(r, R) for the parallel- and antiparallel-spin states of an electron pair. 
Here, r = ri — r2 and R = (ri + r2)/2 are the relative coordinate and the center-of-mass 
coordinate of the pair, the conjugate momenta being k = (ki — k2)/2 and K. The spin- 
resolved pair distribution functions g^jf^i (r) can then be obtained as 



where ((•••)6»)p(fc) represents averages over the scattering angle 6 and over the probability 
p{k) of finding two electrons with relative momentum k. Equation can also be written 
as a double sum of the radial components of the pair wave functions over angular momenta 
and relative momentum, the weight of each term in the latter sum being determined by the 
function p{k) j^. 

Carrying out the angular-momentum expansions of the pair wave functions and perform- 
ing the angular averages leads to the results 



w(r) = ((|vl/^rK(r,R)P)eW) 



(1) 




(2) 



e=l (odd) 



and 
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where Bf^ = 2i + 1 and B^''' = 2. In Eq. Q the sum over £ runs over odd integers, because 
the spinor associated with the IT state is symmetric and hence the coordinate part of the 
wave function is antisymmetric under exchange. The spinor associated with the tJ, state 
has instead no definite symmetry and hence the sum in Eq. ^ runs over all integer values 
of i. The Schrodinger equation obeyed by the radial wave functions ^f^{r) is 



m dr 



m 

iii + 1) (i 



(4) 



0,1,2,.. 



and L2 



2 in 



where Ad = + {D - 1){D - 3)/4, with L| 

{£ = 0, ±1,±2, the negative values of i are accounted for by the choice B^^^ 
Eqs. (0) and Q). In solving Eq. (jlj) we impose as a boundary condition that ^i^{r) tends 
asymptotically to the free solution corresponding to V{r) = 0, except for a phase shift. 

In our calculations the scattering potential V{r) is self-consistently constructed in the 
Hartree approximation. That is, 

V{r) = v{r) +n j d^v'v{\v - r'|) [g{r') - 1] (5) 

where v{r) = e^/r. The system of equations (0)-© is to be solved self-consistently after 
determining the momentum distribution p{k) as a function of temperature, as we do in the 
next Section within a free-electron model. 

The theory satisfies a number of exact properties of g{r), as discussed in more detail in 
Ref. 3: (i) the form of Eq. Q ensures that the relation g-\^{0) = is satisfied, since all 
functions ^i^{r) vanish at the origin for i ^ 0; (ii) the cusp condition on g^i{r) at the 
origin is satisfied if V{r) tends to the bare Coulomb potential for r 0, as imposed in 
Eq. (0); (iii) the normalization of the wave functions and of the probability function p{k) 
ensures that the charge neutrality condition 

n J d^r [g{r) -!] = -! (6) 

is satisfied within a numerical accuracy which is determined by the number of angular- 
momentum states that are taken into account; and finally (iv) the Hartree- Fock expression 



9. 



HF, 



1-5. 



nFY){q) exp (zq ■ r) 



(7) 



is correctly recovered in the weak-coupling limit by setting V{r) = and by describing the 
occupation number of a state with momentum q by the Fermi-Dirac distribution function. 



^FD(g) 



exp 



ksT 



1 



(8) 
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Here, = /i^q^/(2m) is the single-particle energy and hd is the chemical potential, to be 
determined from the mean density n through the relation 



n 



r-\-oo 

/ q''-^ uMq) dq . (9) 
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III. DISTRIBUTION OF RELATIVE MOMENTUM 

The probability distribution p{k) of finding two electrons with relative momentum k is 
given by 

P(^) = ^2 ^2^^D J (^^FD(g)riFD(|q + 2k|). (10) 

The prefactor in Eq. ()10j) is obtained from the normalization condition J^°° p{k)dk = 1, and 
the factor 2k in the argument of the second Fermi-Dirac function takes care of the fact that 
k is the relative momentum of the two colliding electrons. 

At T = the chemical potential is equal to the Fermi energy ep = h'^kp/{2m), with 
kp = (aDrsttB)'^ being the Fermi wave number (as = (97r/4)~^/^ and a2 = 2~^/^), and 
simple analytic expressions are available for p{k) (see Ref. lD and references given there). 
In the limit of high temperature fio approaches the value ksT In [n\^(T)/2] where A(T) = 
^J2tiW' I (mkBT) is the thermal de Broglie wavelength, and Eq. (fTn|l yields the asymptotic 
result 

- ^ '^"-''"m -p (-^) • (11) 

Numerical calculation is necessary at any finite temperature. In D = 2 the specific form of 
the density of single-particle states yields the analytical result 

H2=eF + keT In [1 - exp {-ep/kBT)] . (12) 

The dependence of /id on the reduced temperature 9 is shown in Fig. 1 at = 5 over a 
wide temperature range (recall that room temperature corresponds to O = 0.025). In Fig. 2 
we show the momentum distribution p{k) as a function of k/kp at = 10 in D = 3 and at 
Ts = 5 in D = 2. The effects of temperature are evident from the Figure: (i) p{k) acquires a 
high-momentum tail extending beyond the zero-temperature cut-off at k = kp; and (ii) its 
peak moves to higher momenta from its zero-temperature position at k = kp /2 in D = 3 
and k ^ 0.42A;f in D = 2. 



5 



IV. TEMPERATURE DEPENDENCE OF PAIR DISTRIBUTION FUNCTION 



We report in Figures 3 and 4 our numerical results for the temperature shifts of the 
pair distribution function in the 3D electron gas at = 10 and in the 2D gas at = 5. 
These shifts are qualitatively similar in both dimensionalities and can be said to be relatively 
small, if we bear in mind the magnitude of the unit of measure that we have adopted for 
the reduced temperature. 

Evidently, these results reflect temperature-induced changes in the state of short-range 
order in the gas. Three ranges of values for the reduced distance rkp need to be separately 
examined: (i) the region of the main peak at r^, lying near Tpkp 3.8 in D = 3 and near 
Tpkp «i 3 in = 2, where thermal agitation is seen to reduce the magnitude of two-body 
correlations; (ii) the intermediate region around Tpkp / 2, where the conditioned probability 
of finding a second electron at distance r from an electron at the origin is instead enhanced; 
and (iii) the region of close contact at distances approaching Tpkp ^ 0, where again thermal 
agitation reduces the magnitude of g{r). 

The behaviour of the pair correlations in the first two regions above clearly reflects the 
increase in the mean kinetic energy with temperature, which is associated with the broad- 
ening of the distribution of relative momentum displayed in Figure 2. This enhances the 
tunnelling probability of a second electron into the potential barrier created by the electron 
at the origin and tends to smear out the first-neighbour shell at a given value of the coupling 
strength. Such a reduction in the state of short-range order leads us to expect that the phase 
boundary between the gas and the Wigner crystal will shift with increasing temperature to 
higher values of [i.e. to stronger effective Coulomb coupling). 

The behaviour of qir) in the third region above was first noticed in the work of Macke and 

n 

coworkers relating to the 3D gas in the weak-coupling regime. Their interpretation was 
that the width of the potential barrier surrounding the electron at the origin is also increasing 
with temperature, thus inhibiting deep penetration by a second electron. Our results in Figs. 
3 and 4 confirm this property of the electron gas as persisting in the intermediate-coupling 
regime and in lowered dimensionality. 
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V. SUMMARY AND FUTURE PERSPECTIVES 



In summary, we have solved the two-electron scattering problem in the electron gas at 
finite temperature within the self-consistent Hartree approximation, to display the changes 
in the state of short-range order that are brought about by thermal excitation. 

The results achieved in previous work on pair correlations in the electron gas at zero 
temperature 0, immediately suggest what would be important improvements in the 
present approach, especially in regard to the 2D case and with a view to using the theory 
to construct reliable phase diagrams in a thermodynamic space spanned by density, tem- 
perature, and spontaneous spin polarization. Two directions seem worth exploring: (i) the 
inclusion of a spin dependence of the effective scattering potential, as developed in the work 
of Kukkonen and Overhauser and (ii) the derivation and study of differential equations 
obeyed directly by the pair distribution functions at finite temperature, rather than by their 
angular-momentum components. 

As a final comment, we may recall that in the classical limit the density and temperature 
dependence of g{r) are determined by three- and four-body distribution functions 13|. The 
expression for the temperature derivative of g{r) can then be shown to relate to the heat 
capacity of the classical fluid via a fluctuation formula. The thermodynamics of a Fermi gas 
may also be worth studying from this viewpoint. 
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FIG. 1: Temperature dependence of the chemical potential of a Fermi gas at mean density corre- 
sponding to Ts = 5 in = 3 and D = 2: /Lt/) is plotted in units ep as a function of reduced 
temperature 6. The inset shows an enlargement of the results for 6 < 0.8. 
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FIG. 2: Distribution p{k) of relative momentum k for an electron pair in a Fermi gas at the 
indicated values of the reduced temperature 6, as a function of fc/ZcF for rg = 10 in = 3 (left 
panel) and for = 5 in Z) = 2 (right panel). 
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FIG. 3: Temperature dependence of the pair distribution function in the 3D electron gas at = 10: 
[g{r, @) —g{r, 0)]/© as a function of reduced distance rkp at two values of the reduced temperature 
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FIG. 4: Same as in Figure 3, for the 2D electron gas at rg = 5. 
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